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Abstract 

We establish that the boundaries of the so-called Rost’s reversed barrier are the 
unique couple of left-continuous monotonic functions solving a suitable system of 
nonlinear integral equations of Volterra type. Our result holds for atom-less target 
distributions fi of the related Skorokhod embedding problem. 

The integral equations we obtain here generalise the ones often arising in optimal 
stopping literature and our proof of the uniqueness of the solution goes beyond the 
existing results in the held. 
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1 Introduction 

The celebrated Skorokhod embedding problem formnlated in the 1960’s by Skorokhod |33] 
can be stated as follows: given a probability law p, hnd a stopping time r of a standard 
Brownian motion W snch that Wr is distribnted according to fx. Several solntions to this 
problem have been fonnd over the past fifty years in a variety of settings (see for instance 
Ea, 1271, EH, IB21,1311, HI and [8] among many others; see also the snrvey |26]) and very 
often they hnd important applications ranging from general martingale theory (see for 
instance [J] and [JU]) to mathematical hnance (see [12] or HTl amongst many others, and 
[T8] for a snrvey). 

Here we consider a solntion proposed by Rost ED which can be given in terms of 
hrst hitting times of the time-space process {t, Wt)t>o to a set, depending on /r, nsnally 
called reversed harrier [6] . The bonndaries of this set can be expressed as two monotonic 
fnnctions of time (see i) and we denote them t h->• s±{t) (here we nnderstand s±{t) = 
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it] /^) depending on fi but prefer to adopt a simpler notation as no confusion should 
arise). 

Given a atom-less target distribution ja (possibly singular like the Cantor distribution) 
we use purely probabilistic methods to characterise the couple (s+,s_) as the unique 
solution of suitable nonlinear integral equations of Volterra type. To achieve that we use 
a known connection between Rost’s solution of the embedding problem and an optimal 
stopping problem (see for instance HDl, 1231 and m) so that our contribution becomes 
twofold: on the one hand we provide a new way of characterising and evaluating s± 
and on the other hand we extend in several directions the existing literature concerning 
uniqueness of solutions to certain integral equations. 

To the best of our knowledge until now there exist essentially two ways of computing 
Rost’s reversed barriers: a constructive way and a PDE way. The former approach was 
proposed by Cox and Peskir [S] who approximate the boundaries via a sequence of 
piecewise constant boundaries (in our notation (s!^)n>o)- They prove suitable convergence 
of to as n —)■ oo and for each n they also provide equations that can be solved 
numerically to hnd s!( (however no equations could be provided for s± in the limit). 
In practice for n large enough one can compute an accurate approximation of s±. The 
PDE approach instead relies upon the aforementioned connection to optimal stopping 
as proved in [10] through viscosity solution of variational inequalities. The idea is that 
one should numerically solve a variational inequality to obtain the value function of the 
stopping problem and then use it to compute the stopping set. The boundaries of such 
stopping set indeed coincide with up to a suitable time reversal. An approach of this 
kind was used in Section 4 of CHI, although in the different context of Root’s solutions to 
Skorokhod embedding. 

Here we do use the connection to optimal stopping but instead of employing PDE 
theory we only rely upon stochastic calculus. Our work offers a point of view different to 
the one of |9| and [10] since we directly characterise the boundaries of Rost’s barriers 
as the unique solutions to specihc equations. We do not develop a constructive approach 
nor we use an indirect approach via variational inequalities. The computational effort 
required to evaluate our equations numerically is rather modest and hinges on a simple 
algorithm which is widely used in the optimal stopping literature (see Section]^. It is 
worth mentioning that McConnell [23] provided a family of uncountably many integral 
equations for which appear to be very different from ours as they involve at the same 
time the boundaries and (in our language) the value function of the optimal stopping 
problem. Solvability of these equations seems to be left as an open question and we are 
not aware of any further contributions in the direction of Ending its answer. 

To clarify our contribution to the literature on integral equations and optimal stop¬ 
ping we would like to recall that nonlinear equations of Volterra type arise frequently in 
problems of optimal stopping on finite time horizon. Such equations are indeed used to 
compute the related optimal stopping boundaries (see for a collection of examples) 
and their study in this context goes back to McKean and Van Moerbeke [35] among 
others (see also [5], [19], [21] and [25] for the American put exercise boundary). Volterra 
equations are also widely studied in connection with the Stefan problem in PDE theory 
(see for example [1] and Chapter 8 in [14jh 

The particular equations that we consider here (see Theorem 2.3 below) are more 
general than the ones usually addressed in optimal stopping literature and in PDE theory. 
All examples we know from optimal stopping involve integrals with respect to the Lebesgue 
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measure on M+ x M. Our equations instead involve integrals with respect to a product 
measure dt x [v—gi) [dx) on M_|_ x M, where g, and u are rather general probability measures. 
Only when both g and u have density with respect to the Lebesgue measure we can rely 
npon the existing results. 

As for the integral eqnations used in PDE theory in the context of the Stefan problem, 
there is no direct way of employing the related existence and uniqueness results (see for 
example Chapter 8 in and compare to our (2.8)). This is due to the need for a regularity 
of the terms in the equation that cannot be achieved when considering Rost’s solutions 
to the Skorokhod embedding problem. 

In this sense our Erst contribution is to establish existence of solutions for such broader 
class of integral equations. Secondly, the question of uniqueness of the solution was a long 
standing problem for those equations normally addressed in optimal stopping theory. It 
found a positive answer in the paper by Peskir [2S] where uniqueness was proved in the 
class of continnous fnnctions. Here we prove that such uniqueness holds for a more general 
class of integral equations and, more interestingly, we prove that the solution is indeed 
unique in the class of right-continuous fnnctions, thus extending Peskir’s previous results 
(note that onr result obviously holds for all the cases covered by [2H])- We remark that 
although the line of proof is based on a 4-step procedure as in [2S], the greater generality 
of our setting reqnires new arguments to prove each one of the four steps so that our proof 
could not be derived by the original one in |28j . 

Before concluding this introduction we would like to point out that our work is some¬ 
how related to a recent one by Gassiat, Mijatovic and Oberhauser |T5]. There the authors 
consider another solution of the Skorokhod embedding problem, namely Root’s solution 
[32] . which is also given in terms of a barrier set in time-space. The boundary of Root’s 
barrier is expressed as a fnnction of space x h-)■ r{x) rather than a fnnction of time. Un¬ 
der the assnmption of atom-less target measures g (as in our case) [15] prove that the 
boundary r( •) solves a Volterra-type eqnation which is reminiscent of the ones obtained 
here. Uniqueness is provided via arguments based on viscosity theory in the special cases 
of laws g’s which produce continuous boundaries r( ■). This additional requirement ex¬ 
cludes for instance continuous distributions singular with respect to the Lebesgue measure 
(e.g. Cantor distribution) which are instead covered in our paper. 

It is important to notice that although Rost’s and Root’s solutions are conceptually 
related there seems to be no easy way of obtaining one from the other. Therefore even if 
equations found here and those in [T5] may look similar we could not establish any rigorous 
link. Rost’s barriers enjoy nicer regularity properties than Root’s ones (e.g. monotonicity 
and left continuity) and this allowed us to prove uniqueness of solution for the integral 
equations even in the case of continuous singular measures g. For the same reason the 
Volterra equations that we obtain here are relatively easy to compute numerically; indeed 
we can rnn our algorithm forward in time starting from known initial values s±(0) of the 
boundaries at time zero. The equations obtained in [15] are instead more complicated as 
there are no clear boundary conditions on r( ■) that can be used in general to initialise the 
algorithm. In fact the numerical evaluation of the boundaries in ng is restricted to the 
class of g’s that produce symmetric, continuous, monotonic on [0, oo], maps x i—)■ r{x). 
Finally the existence results in our paper and in [TS] are derived in completely different 
ways. 

The rest of the paper is organised as follows. In Section]^ we provide some backgronnd 
material on Rost’s barriers, state our main result in Theorem [2.3 and recall the link to 






Integral equations for Rost’s barriers 


4 


a suitable optimal stopping problem which is needed for the proof. Section is entirely 
devoted to the proof of Theorem |2.3 Finally, in Section we address the numerical 
solution of the integral equations. 


2 Setting, background material and results 

Let P) be a probability space, B := {Bt)t>o a one dimensional standard Brown¬ 

ian motion and denote {fFt)t>o the natural filtration of B augmented with P-null-sets. 
Throughout the paper we will equivalently use the notations E/(i?f) and f ■ 

M —)■ M, Borel-measurable, to refer to expectations under the initial condition Bq = x. 
We also use the notation X for a random variable X with law /i. 

Let fi and z/ be probability measures on M. We denote by F^{x) := /r((—oo,a;]) and 
F^{x) := z/((—cx), x]) the (right-continuous) cumulative distributions functions of fi and 
u. Throughout the paper we will use the following notation: 

a+ := sup{x G M : X G supp z/} and a_ := — inf{x G M : x G suppz/} (2.1) 

/i+ := sup{x G M : X G supp/i} and := — inf{x G M : x G supp/x} (2.2) 

and for the sake of simplicity but with no loss of generality we will assume a± > 0. We 
also make the following assumptions which are standard in the context of Rost’s solutions 
to the Skorokhod embedding problem (see for example [9] and in particular Remark 2 
therein). 

(D.l) There exist numbers 6_|_ > a+ and b_ > a_ such that (—&_, &+) is the largest interval 
containing (—a_,a+) with /i((—6_,6+)) = 0; 

(D.2) If 6+ = a+ (resp. b_ = a_) then /x({6_|_}) = 0 (resp. /x({—6_}) = 0). 

It should be noticed in particular that in the canonical example of v{dx) = 6o{x)dx the 
above conditions hold for any fr such that /x({0}) = 0. We also observe that assumption 
(D.2) is made in order to avoid solutions of the Skorokhod embedding problem involving 
randomised stopping times whereas assumption (D.l) guarantees that at each time t > 0 
the t-section of the Rost barrier is the union of at most two connected components. 

Now we recall a particular result recently obtained by Cox and Peskir j9] in the context 
of a more general work on Rost’s solutions of the Skorokhod embedding problem. Here we 
give a different statement for notational convenience and the next theorem summarises 
P Thm. 1, Prop. 9 and Thm. 10]. 

Theorem 2.1 (Cox & Peskir |9]). Let W'' := {W^)t>o be a standard Brownian motion 
with initial distribution z/ and let (D.l) and (D.2) hold. There exists a unique couple of 
left continuous, increasing functions s+ : M+ —)■ IR+ U {+C)o} and s_ : M+ —)• IR+ U {+C)o} 
such that if we define 

(T* := inf [t > 0 : Wf ^ ( — s_(f), s+(f)) } (2.3) 

then ~ /X. 

Moreover, letting As±{t) := s±{t+) — s±{t) > 0, for any t >0 such that s±{t) < -foo 
it also holds 


As+(f) > 0 /x((s+(f),s+(f+))) = 0 

As-{t) >0 44. t^{{~ 5-(i+), — 0- 


(2.4) 

(2.5) 
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We point out that the above theorem is stated in [9l Thm. 1] for the special case 
v{dx) = 6Q{x)dx. Remark 2 in [S] clarihes that the result and its proof remain identical 
for any initial distribution satisfying (D.l) and (D.2). Alternatively the same result may 
be found in m Thm. 2.3] where a different proof relying upon optimal stopping theory 
is provided for general u. 

Since are the boundaries of the Rost barrier associated to /i, it is clear that flat 
portions of or s_ must correspond to atoms of /i. Hence the next corollary follows. 

Corollary 2.2. If the measure p is atom-less then the boundaries s+ and s_ are strictly 
monotone. 


In preparation for the statement of our main result let us denote 

p{t,x,s,y) Q- 2 (s-t) ^ for f < s, X, 1 /e M (2.6) 

^27r(s-f) 

the Brownian motion transition density and let (L^)s>o be the local time of R at a point 
1/ G M. 

As in Section 4 of m we introduce the (generalised) inverse of dehned by 

{ inf{f > 0 : —S-{t) < x}, x < —s_(0), 

0, X e (-s_(0),s+(0)), (2.7) 

inf{f > 0 : s+(f) > x}, x > s+(0). 

Note that for f > 0, one has x G (—s_(t), s+(f)) if and only if (p(x) < t. Moreover (p is 
positive, decreasing left-continuous on ]R_ and increasing right-continuous on IR+ (hence 
upper semi-continuous). 

Now that we have introduced the necessary background material we can provide the 
main result of this work, i.e. the characterisation of the boundaries as the unique 
solution of suitable nonlinear integral equations of Volterra type. Its proof will be given 
in the next section. 


Theorem 2.3. Assume that y is atom-less. Then for any T > 0 the couple (s+,s_) is 
the unique couple of left-continuous, increasing, positive functions that solve the system 



i's+{T-u) 


'-s-{T-u) 


p{t, ±s±(T -t),u, y) (z/ - /i) {dy) du = 0, 


f G [0,T) 


( 2 . 8 ) 


with initial conditions 5±(0) = b±. Equivalently we may express (2.8) in terms of cp as 


(f- n)(dy) = 0 , X 


(2,9) 


Remark 2.4. It is important to notice that (2.8) reduces to a single equation if supp{y}n 
M_ = 0 (hence S-(-) = +C)oj or supp{y} fl M+ = 0 (hence s+(-) = +C)oj. Figures 2 and[^ 
illustrate numerical solutions of (2.8) also in cases with a single boundary. 
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Since the proof of the above theorem relies on the analysis of a problem of optimal 
stopping, then we introduce such problem here before closing this section. For 0 < T < 
+00 and (t,x) G [0,T] x M we denote 


G{x)-.=2[ {F,{z) - F^{z))dz 
Jo 

and dehne 


( 2 . 10 ) 


V'^{t,x) := sup Ea;G{Br) 

0<T<T-t 


( 2 . 11 ) 


where the supremum is taken over all (J^i)-stopping times in [0,T — t]. Here the upper 
index T in the dehnition of the value function is used to account for the time horizon 
explicitly. As usual the continuation set Ct and the stopping set Vt of (2.11) are given 
by 


Ct := G [0,T] x M : V^{t,x) > (2.12) 

Vt := {{t,x) G [0,T] X M : V^{t,x) = G'(x)} (2.13) 


and clearly 


Vt — ([0, F] X M.) \ Ct- 


Throughout the paper we will often use the following notation: for a set A C [0,T] x M 
we denote Af] {t < T} := {(t, x) G A : t < T}. 

Problem ( |2.11[ ) was addressed in m where it was shown that the geometry of the 
continuation and stopping sets can be characterised in terms of the boundaries of the 
reversed Rost’s barrier. Here we summarise those results in the following theorem under 
assumptions (D.l) and (D.2). 

Theorem 2.5 (De Angelis |II]). For each T > 0 it holds 


Ct = {(^,a;) G [0,T) X M : a; G ( — s_(T — t),s+(T — t))} 


and the smallest optimal stopping time of problem (2.11) is given by 
Tj'{t,x) = inf {s G [0,T - f] : (f + s,x + Bg) ^ Ct} 


(2.14) 

(2.15) 


for {t, x) G [0, T] X M. Moreover the functions fulfill the following additional properties: 
{T — f) > h± for f G [0, T) and (0+) = h±. 


For T > 0 it is convenient to introduce the function 


U'^it, x) := V'^(t^ x) — G{x), (t, a;) G [0, +oo) x 


(2.16) 


Applying Ito-Tanaka-Meyer formula to G inside the expectation of (2.11) we hnd an 
alternative representation for U'^ which will also be useful in the rest of the paper, i.e. 

U'^{t,x)= sup [ - fi){dy). (2.17) 

0<T<T-t 7r 


In order to prove that the couple (s+, s_) solves the integral equations (2.8) we use a 
characterisation of the function U'^ which is of independent interest. We provide this in 
a separate proposition whose proof will be given in the next section. 
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Proposition 2.6. Assume that fi is atom-less. Then for any T > 0 the following equiv¬ 


alent representations of (2.16) hold 


U'^it, x) = 


rs+{T-u) 

/ P{t,x,u,y){u - fi){dy)ds 

J-s.{T-u) 

(2.18) 

{y ■ ‘p(y)<T-t}^xLT-t-ip{y) ~ 

(2.19) 


for {t,x) e [0,T] X 


3 Proof of the main theorem 


In this section we provide the proof of Theorem 2.3 This requires to show hrst that the 


value function V'^ of (2.11) is in the whole space and then to introduce a free boundary 
problem. The latter allows us to obtain the representation formulae of Proposition 2.6| 
which hnally lead to the result in Theorem 2.3 These steps will be taken in three separate 
sub-sections. 


3.1 regularity of the value function and a free boundary problem 

We start by recalling that Proposition 3.1 in m guarantees that 

e C'([0,T] X M). (3.1) 

Moreover it has been shown in m , and it is one of the paper’s main technical results, that 
the time-derivative Vf’" of the value function is indeed continuous on [0, T) x M. However 
it is also mentioned there that when p has atoms the continuity of Vj'{t, ■) could break 
down across the boundaries (let alone continuity of ■)). For this reason here we 
focus only on continuous target measures p. 

First we summarise the results of Lemma 3.7 of m which will be needed in the 
following Proposition. 

Lemma 3.1. For any (t,x) G dCx fl {t < T} and any sequence (tn,Xn)n C Ct such that 
(tn,Xn) —t (t,x) as n ^ oo it holds 

Mm rj'itn, Xn) = 0, P — a.s. (3.2) 

n^O 

Next we provide the regularity of V'^. 

Proposition 3.2. Assume that /i is atom-less. Then for any T > 0 it holds V'^ G 
C'^([0,T) X M). 

Proof, (a). As already mentioned one can check in] for V/ G C([0,T) X M). 

(b). Here we prove that G ^([OjT) x M). The claim is obvious in the interior of 
the continuation set Ct where it is well known that V'^ G (see [2^ Sec. 7.1]). In the 
interior of the stopping set Vt one has = G and Vj = G' = 2{F^ — F^). The latter is 
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continuous in Vt fl {t < T} because for all t G [0, T) one has s±{t) > b± due to Corollary 


2.2, and therefore 


Fu{x) = 0 for a; < —s_(f) and Fi^{x) = 1 for a; > s+(t), 

for all t G [0,T) such that s+(t) < +cxd and/or s_(f) < +cxd. 

Then we must only consider points along the two boundaries. Let {t, x) G OCt fl {f < 
T} be arbitrary but hxed and let {tn,Xn)n C be a sequence of points converging to 
(t,x) as n — >■ cxD. Fix an arbitrary n G N and for simplicity denote := rJ{tn,Xn) (see 

Since stopping at is optimal for the initial conditions (tn,Xn), then for h > 0 such 
that (tn, Xn + h) E Ct we use the mean value theorem to obtain 


XnFh)- Xn) ^ 1 


h 




G(X^-+f‘) - G(XZ‘) = EG'(Jf) 


f ( ^h\ 


(3,3) 


with G (X/", P-a.s. On the other hand by the same arguments we also obtain 


V'^{tn,Xn) -V'^{tn,Xn- h) 

h - h 


G(X/:) - G(X/:-") = EG\^^) 


llGh\ 


(3.4) 


with ^2 ^ P-a-s. 

From continuity of the sample paths we get that and converge to X/" 

as h —)■ 0, P-a.s., hence ^2 as h —)• 0 as well. Recall that I 4 exists inside Ct 

and G is Lipschitz on M, so that dominated convergence and (3.3) and (3.4) allow us to 
conclude 


Vj{tn,Xn) = EG'{X::). 


(3.5) 


Now we use continuity of sample paths and Lemma 3.1 to obtain X/" —>■ x as n —)■ cx). 
Hence another application of dominated convergence and (3.5) give 

lim Vj{tn,Xn) = G'{x). (3.6) 

n^oo 

Since the point (f,x) and the sequence (tn,Xn) were arbitrary we conclude that is 
continuous across the boundaries and hence everywhere in [0, T) x M. □ 


It follows from the standard theory (see for instance [29l Sec. 7.1]) that (3.1) and the 
strong Markov property imply that V'^ G in Ct and it solves 


(rt + ivS)(i,i) = o, 

V'^{t, x) = G{x), 


{t, x) G Ct 
{ t,x) G OCt 


(3.7) 

(3.8) 


where Ct may be expressed in terms of the curves s±{T — ■) as in (2.14). It now follows 
from Proposition 3.2 that if fi is atom-less then is also globally G^ in [0,T) x M and 
two more boundary conditions hold as well, i.e. 


= G"(x), 

yt^{t,x) = 0, 


(f, x) G OCt n {f < t} 
{t, x) G OCt n {f < t}. 


(3.9) 

(3.10) 
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It is important to notice that is continuous on [0,T) x M even without assuming 
that fi is atom-less (see HD Prop. 3.15]). Then is continuous on the closure of the 


set Ct fi {t < T — e} for any e > 0, due to (3.7) and continuity of V^. In Vt we have 
yT _ qt therefore 14^(t, dx) = G'\dx) exists as a measure therein. 

The above considerations imply that the function U'^ = — G (see (2.16)) is con¬ 


tinuously differentiable in t and ■) dehnes a measure on M for each hxed t G [0, T). 


Notice that if // is atom-less then thanks to (3.9) such measure charges no mass at 
dCx n {f < T}, hence the above boundary value problem can be stated for U'^ as 


{u[ + = 

U^{t, x) = 0, 
Uj{t,x) = Uj{t,x) 


(z/ — ^){dx), 


= 0 , 


{t, x) G Ct 
( t, x) G OCt 
{ t, x) G OCt n {f < T} 


(3.11) 

(3.12) 

(3.13) 


where (3.11) holds in the sense of distributions. 


3.2 A representation of the value function 


Here we prove Proposition 2.6 The integral representation (2.18) for easily yields 


the integral equations for the boundaries and it is therefore instrumental for our main 


result, i.e. Theorem 2^ Moreover some of the concepts introduced in the proof below 
will be needed when proving the latter theorem. 


Proof of Proposition 2.f^ For the sake of this proof it is convenient to denote b±{t) : — 
s±{T — t) for t G [0, T] and use h± instead of s±. Similarly we drop the apex T and denote 
U = U^. 

The representation (2.18) for U'^ is obtained by an application of an extension of Ito- 


Tanaka-Meyer formula to the time-space framework which is inspired by [T]. However 
since the result of [H Prop. 4] cannot be applied directly to our setting we provide a full 
proof for completeness. 

We start by recalling that for each t G [0, T) the map x h-)■ Uxit, x) may be seen as the 
cumulative distribution function of a signed measure on M with support on [—&_(f), b+(t)]. 
In fact we set U^xit, {a,b]) := Ux(t,b) — Ux(t,a) for any interval (a, 6) C M (recall that 


Fy is only right continuous) and we notice that due to (3.13) Uxx{t, {±6±(t)}) = 0 for all 
t G [0,T), i.e. Uxx{t,dx) has no atoms at the boundary points b±{t). It then follows 

Uxx{t,dx) = l{x&(-b_{t),b+{t))}(Vxx{t,x)dx -2{iy - n){dx)^, {t,x) G [0,T) x M. (3.14) 

We continuously extend U to by taking U(t,x) = U{0,x) for f < 0 and U(t,x) = 
U (T,x) = 0 for f > T. Then we pick positive functions f,g^ G"^ (M) such that 


f{x)dx = / g{x)dx = 1 


and for u G N we set fn{,x) := nf{nx), gn{x) = ng{nx). We use standard mollihcation to 
obtain a double indexed sequence of functions (f7n,m)n,meN2 C C'“(M^) dehned by 

Un,m{t,x)-.= f f U{s,y)fn{x - y)gm{t - s) dyds, {t,x) eMf. (3.15) 
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To avoid technicalities related to the endpoints of the time domain £x [t, x) G (e, T—e)xM., 
£ > 0 arbitrary, then by applying classical Dynkin’s formula for each n and m we hnd the 
expression 

r 1 

= + 1 ( 5 , )dsj. (3.16) 

By the standard properties of mollihers it is not hard to see that taking limits as n, m —)■ oo 
and as e —)• 0 one has Un,m{t, x) —)■ f/(f, x), E [Un,m{T — e, E[U (T, B!fl’)] = 0. 

It only remains to analyse the integral term. Note that when extending U to we 
have implicitly taken b±{t) = &±(0) for t < 0 and b±{t) = 0 for t > T. Note as well that 
there is no loss of generality assuming that the index m in (3.16) is large enough to have 
Qm supported in {—e, +e:). Then for (s, z) E (t,T — e) x R. it holds 

r-fe+(l)) 


i£,U^,m){s,z) = 
'^Un,m) {s 1 z) = 



u(v, y)fr, (z - y)gm{s - v) dy dv 



-6- (d) 

fe+h) 

-b- (d) 


U(v,y)f„(z - y)g'„(s - v) dydv 


(3.17) 

(3.18) 


where the integrals in dv must only be taken for v E {s — e,s + s) C (0, T). 
Use integration by parts twice, the smooth-fit (3.13), and (|3.14[) to obtain 


{£,Un,m){s,z) = 


pb+iv) 

' —b- (v) 
rb+{v) 

' —b- (v) 


fn{z - y)U^xiv, dy) jgm(s - v) dv 


(3.19) 


fn{z - y){Vxx{v,y)dy - 2(n - y){dy)) )gm{s - v) dv 


and similarly 


{dt^n,m} (s, z) — 



b+{v) 


-b-{v) 


Ut{v, y)fn{z - y)gm{s - v) dydv. 


(3.20) 


Since Ut = Vt and (3.7) holds, then adding up (3.19) and (3.20) (where now 2 ; = 5*’"^) the 
above gives 


r-T-e 


Jt 

= - E 


ilU^,m + li^Un,m){s,Bl’^)ds 

rb+iv) 


(3.21) 



pT-e 


T-e 


'-b- H 
rb+H 


fn{Bf -y){iy- y){dy) jgm(s- v) dv 


'-b-(v) 


fn{z - y){v - y){dy)^gm{s - v) dn|p(f,a;, s,z)ds dz. 


Recall that supp^f^ C (—£, -1-e). For (s, z) E {t,T — e) x R we set 


An{v,z) := 


i>b+{v) 


fn{z - y){iz - fi){dy), vE{s-e,s + e) 


(3.22) 


J-b-{v) 

An,m{s,z)\= / kn{v,z)gm{s -V)dv 


(3.23) 
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so that we rewrite (3.21) in the form 

rT-e 

1 





T-e 


An,m(s, z)p{t, X, S, z)ds dz. 

(3.24) 


For hxed n and for any G M and s E {t,T — e) the map v h-)■ A„(n, z) is bounded and 
continuous on {s — e,s + e). This follows since fi is flat at possible jumps of h± (cf. (2.4), 
(2.5) of Theorem 2.1 and recall that h±{-) = s±{T — ■)) and by observing that 



y)T^{.dy) 



y)T^{.dy). 


Notice also that gm{s — ■) converges weakly to (5(s — •) as m —)• oo and hence for all z G M 
and s G (t,T — e) we get 


lim An,mis,z) = An{s,z). 


(3.25) 


For hxed n it holds ||/n||oo < Kn for suitably large Kn > 0 and hence it also holds 
I An,m(s, z) I < 2Kn Uniformly for (s, z) G {t,T — e) x M. We keep n hxed and take limits 
as m —)■ oo, then dominated convergence and (3.25) give 


lim 

m^oo 



T-e 


An,m{s, z)p{t, X, S, z)ds dz = 



T-e 


An{s,z)p(t,x,s,z)dsdz. (3.26) 


In order to take limits as n —)■ cx) we analyse the the right-hand side of the expression 
3 ve. In particular we recall from (2.7) that 

y G (-6-(s), 6+(s)) -^y E (-S-(T - s), s+(T - s)) ip{y) <T -s. (3.27) 


Hence Fubini’s theorem and the occupation formula for local time give 

rT-e 

A„(s, z)p(t, X, s, z)ds dz 



=E 


t.X 


rT-e rh+is) 


^ J t J —b- (s) 

rT-e 


fn{Bs - y){iz - p){dy) ds 


= E 


t^X 


t{s<T-^(y)}fn{.Bs - y)ds {u - p){dy) 


= E 


t,X 


gT-e)A{[T-^iy)]Vt) 


fn{Bs-y)ds {iz - p){dy) 


Jt 


fn{z - y)Et,a: [Lljr_,)^(^[T-y,^y)]^t)]dz) {iZ - p){dy) 


(3.28) 


where we recall that {Ll)s>o is the local time of the Brownian motion at a point x G M. 

With no loss of generality and thanks to the dehnition of /„, we may assume that 
there exists a compact iF C M such that supp /„ C 7F for all n. For hxed t and x the 
quantity \LlT-e)A([T-y:>{y)]\/f^ bounded uniformly with respect to all y, x G M such 
that y — z E K (use for example Ito-Tanaka’s formula to verify the claim). Moreover, 
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for fixed t, x and y the map i—)■ 3 , is bounded and continuous for 

z G M such that z — y G K. Finally, for each y E M. one has fn{- — y) —t (5(- — ?/) weakly 
as n —)■ cx). 

The above considerations allow us to use dominated convergence (with respect to the 
measure v — y) and (3.28) to obtain 


lim 

n^oo 



T-e 


An{s,z)p{t,x,s,z)dsdz = / Et,x[L\T-e)Ai[T-^iy)]vt)]i’^ - (3-29) 


Now, we collect (3.16) and (3.24), we take limits as m —)■ cxd hrst, then as n —)■ cxd and 
hnally as £ —)■ 0. Then using (3.26) and (3.29) we obtain 


U{t,x)= I - y){dy) 


(3.30) 


/ ^{y-‘piy)<T-t}^x[L^_t_y,ty\\{^ A){dy). 


The equivalent expression (|2.18|) can be obtained by recalling the well known representa- 

(3.31) 


tion for the expectation of the Brownian local time 

rT-ip{y) 


p{t,x,s,y)ds. 


Then by using once again (3.27) we obtain 


U{t,x)= l^s<T-^iy)}P{t, X, s, y)ds^ (z/ - /i) (dy) 


(3.32) 


l{ye(-6-(.).fe+(.))}P(i, X, s, y)ds ) (z/ - p){dy) 


and hence ( |2.18[ ) follows by a simple application of Fubini’s theorem and by putting 
b±it) = s±{T-t). □ 


3.3 Uniqueness of the solution for the integral equations 

We can now prove Theorem 2.3[ 


Proof of Theorem 2^. We develop the proof considering a setting with two boundaries, 
i.e. /i supported on both M_ and IR+. A simpler version of the same arguments holds in 
the cases with a single boundary. 


First we notice that the system (2.8) of integral equations is easily obtained by taking 
X = ±s±(T — t), t < T in (2.18) and observing that U^(t, ±§±(7" — t)) = 0. For the rest 
of this proof it is convenient to denote b±(t) := s±{T — t) fort G [0, T] and use b± instead 
of s±. Also we drop the apex T and denote U = U'^. 

The proof follows a contradiction argument. Let us assume that there exists a couple 
of functions (r_,r+) with the following properties: r± : M+ —)■ M_|_ U {+cx)}, increasing. 


left-continuous and such that (r_,r+) solves ( 2 . 8 ) with terminal conditions r±( 0 -|-) = b 


->±- 


Notice that by monotonicity r±{T — t) > a± for t G [0,T). As in (2.7) we can define a 
generalised inverse of r± and we denote it by V’ : 1 ^ —t [ 0 , cx)). 
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Again it is convenient to consider a time reversal and denote c±{t) := r±{T — t), 


t G [0,T]. Motivated by (2.19) let ns define a fnnction 


U%t,x) := j {t,x) G [0,T] x M (3.33) 

J M 

and notice that for we simply have If’ = U (observe that ^{y:ip{y)<T-t}LT-t-‘)b(v) ~ 


^?T-t-V>(y))+ 


). By assnmption (c_,c+) solves ( |2.8 ) and hence 

U'’{t,±c±{t)) = 0, tG[0,T), and U'’{T,x) = 0, xE 


(3.34) 


Using that = E^\B(T-t-ij{y))+ - y\ - \x - y\ it is not difficnlt to prove 

that G U([0,T] X M). Now we want to show that the process Y” defined as 


Yr.= U'’{t + s,Bs)+ I sE[0,T-t] 


is a Paj-martingale for any a; G M and t G [0,T). By the definition (3.33) it follows 


ExU {t S, Bg) — Ex / l.[t+s<T-^(y)}EBs[L^(^T-(t+s)-ip(y))\(^ h')('^l/) 

Jm. 

Strong Markov property and Ito-Tanaka’s formnla give 


(3.35) 


(3.36) 


Eb LH = Ea 


Bs-\-u y 




\Bg -y\ = Ex{Ll^J\Bs) -LI- - a.s. (3.37) 


for all M > 0 and with {Ms)s>q a P^^-martingale with Mq = 0, P^j-a.s. In particnlar for 


M = T — (f + s) — V’( 2 /), by Fnbini’s theorem and (3.37) we obtain 

ExU%t + s, Bg) = f l{g<T-t-^{y)}Ex [h^T-t-p(y)) - - y)idy). (3.38) 

M 

Now we add and snbtract U^(t,x) from the right-hand side of the above expression, use 


(3.33) and simplify the indicator variables to get 


ExU^it + s,Bg) =U%t,x) 

J M 

~ / ^{s<T-t-'ilj{y)}Ex Ll{v- y){dy) 

JM. 

=U {t,x) — I l.^Q^T-t-‘)p{y)<s}ExLjrp_^_^^y^^^{i' — y)(dy) 


(3.39) 


/ l{s<T—4—p(j/)}Ea;h^(i/ y){dy^ 

JR 

=U^{Lx) - Ex - y){dy). 

Therefore the process Y^ is a continuous martingale as claimed and from now on we will 
refer to this property as to the martingale property of Lf^. Further, we notice that in the 
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particular case of c± = b± we obtain what we will refer to as the martingale property of 
U. Then for any {t, x) G [0, T] x M and any stopping time r G [0, T — f] it holds 


U{t,x) = £^ U{t + T,Br)+ / 

U^{t,x) = E^ f/"(t + r,5^)+ f - p,){dy) 


(3.40) 

(3.41) 


The equations considered by 
Peskir can be seen as a special case of (2.8) and (2.9) and here we obtain uniqueness in a 


We now proceed in four steps inspired by Peskir 


stronger sense, i.e. in the class of right-continuous functions rather than continuous ones 


as in 



Figure 1: Examples of Brownian sample paths prior to the hrst hitting time r to the 
boundary c+. 


Step 1. First we show that 

x) = 0, for t G [0, T), X G (—cxD, — c_(f)] U [c+(f), oo). (3.42) 

We £x f G [0,T) and with no loss of generality it is sufficient to consider x > c+(t) 
since the proof for x < —C-{t) follows from analogous arguments. We set r ;= inf{s G 
(0,T — t] : Bs < c+(t -|- s)} and note that since s h->■ Bs{uj) — c+(t -f s) is positive for 
s < T, right-continuous, with only upwards jumps, then Bs{oj) — c+(t + s) must cross 
zero in a continuous way and it must be B^ = c+(t -f r), P^j-a.s. Then (|3.34) implies 


lB’{t -|- t,Bt-) = 0, Paj-a.s. and it follows from (3.39) that 


f/"(f,x)= / E,L^^(^_^_^(^))+(z/-/i)(di/). 


(3.43) 


The support of y contained in [6_|_, c_|_(t)), P 2 .-a.s. by dehnition of r 

and by observing that 


{z G M : 'ip{z) <T — t} = (—c_(t),c+(f)) 


(3.44) 
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with c+(T) = 6_|_. Since c+ is monotonic decreasing, for any y > 0 it is not hard to see 
that: 


i) if t{u)) <T — t — 4’{y), then Bs{uj) > y for s < t{uj)] 
if) if t{uj) > T — t — i/j{y), then Bs{uj) >yioTS<T — t — t/jiy)] 


hence = 0, Pa;-a.s. and (3.42) follows (see also Figure 1). Since U is non 

negative by definition (see (2.16)) then it also follows 


0 = U^(t, x) < U{t, x), for f G [0, T), a; G (— cxd, —c_(f)] U [c+(f), oo). 


(3.45) 


Step 2. Next we prove that <U on [0,T] x M. Since U^{T,x) = U{T,x), a: G M 
and (3.45) holds it only remains to consider t G [0,T) and x G (— c_(t), c+(t)). Fix any 
such (f, a;) and denote a := inf {s G [0, T — f) : Bg < —C-{t + s) or Bg > c+{t + s)}. Notice 


that due to (3.42) and monotonicity of c+ and c_, we have U‘’{t + a^B^-) = 0, Pa,-a.s. 


Then from (2.17) and the martingale property of we obtain 


U{t,x) -U\t,x)> f y){dy) 

•J M 

>- f E4L> - (,))*] 

J M 


(3.46) 


since u h->• L^{uj) is increasing for each ?/ G M. We recall (3.44) and claim that 


TV — TV 


Pa,-a.s. for y G (-c_(f), c+(f)). 


(3.47) 


The latter is trivial if cr{u;) < T — t — 'ipiv)- Consider instead cr{u) > T — t — 'ip{y) and 
with no loss of generality assume ?/ > 6+. Monotonicity of c+ implies that Bs{uj) < y 
for s G (T — t — 'ip{y),a{u)) (see also Figure 1) so that the local time cannot increase 
after time T — t — 'ip{y) and therefore the claim holds. The case of y < —b_ follows from 
symmetric arguments relative to c_(-). 


From (3.47) and observing that y is supported in (—c_(f), c+(t)), P^j-a.s., we 


conclude that the right hand side of (3.46) equals zero and <U 


Step 3. Now we aim at proving that c±{t) < b±(t) for t G [0,T). Let us assume 
that there exists t G [0,T) such that c+(t) > 5+(t) and take x > c+(f). We denote 
t' := inf{s G [0,T — f) : Bg < b+(t + s)} and from the same arguments as those 
used in step 1 above we have B^-' = b+(t + r') and U(t + t',Bt-i) = 0, Pj,-a.s. Since 
X > c+(t) > 5+(t) we have also U'’{t,x) = U{t,x) = 0 and from step 2 above it follows 


0 = U{t + t\Bt') > U’’(t + r',i?T-')) Paj-a.s. Subtracting (3.40) from (3.41) now gives 


'i^l'A{T-t-'iP(y))+ ^l'AiT-t-y}{y)) + )(^ y){dy) > 0. 


Due to monotonicity of 6+, an argument similar to the one that allowed us to show that 
the right-hand side of (3.43) vanishes, here implies 


= 0, P 3 .-a.s. Therefore 


from the last inequality it follows 


0 < 


F r2/ 

'^x-^^T'A{T-t-^P{y))+ 


{v - fi){dy) = - j < 0, (3.48) 
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where we have also used that y G supp z/ ^ = 0, Pa;-a.s. since x + > a+ 

for s < t' (recall that &+(■) > a+ in [0,T)). 


Now (3.48) implies that 

for /i-a.e. y G M, = 0, P^^-a-s. 


(3.49) 


Since y h->■ is supported in [ 6 +,c+(t)), P 3 ;-a.s. and x > c+(t), then (3.49) 

also implies that 


T A{T - t- ip{y))^ < Cy, Px-a.s., for p-a.e. y G 


(3.50) 


with C^y := inf{s > 0 : = y}- Next we show that (3.50) leads to a contradiction. 

Let / C ( 6 +(t), c+(t)) be an open interval. For any y E I, monotonicity of 6 + implies 
that t' > Cy, 

(T 


a.s. and therefore it follows from (3.50) that for /r-a.e. y E I it must be 
P 3 .-a.s. However, for any hxed y G /, by right-continuity of c+ there 


■t-i^{y))^ < Cy, 

must exist £ > 0 such that c+(t -|- s) > y for s E [0, e:] and equivalently 'ip{y) < T — {t + e). 
For such y and e this leads to {T — t — < Cy ^ Cy > P^-a.s., which is clearly 

impossible since Px{,Cy < ^) > 0. Therefore c+(t) < h+{t) for all t E [0,T). 

To prove that c_(t) < 5-(t), t G [0,T) we assume that there exists t G [0,T) such 
that c_(t) > 5-(t); then we pick x < —c_{t) and follow arguments as above to show a 
contradiction. 

Step 4- At this point we prove c± = b±. Again it is enough to give the full argument 
only for the upper boundaries as the case of the lower ones is similar. 

We assume that there exists t G [0,T) such that c+(t) < b+{t), then by dehnition 
of 6+ there must exists 5 > 0 and x E (c+(t), 6+(t)) such that U{t,x) > 5. We denote 
a' := inf{s G [0,T — t) : Bg < —b_{t + s) or Bg > 5+(t -|- s)} and we recall that a' is 


optimal for U(t,x) (i.e. a' = r* as in (2.15)), so that (2.17) gives 


U(t,x)^ / - ii)(dy). 


(3.61) 


Since c± < b± (from step 3 above) and U^{T,x) = 0, it follows that U'’(t,x) = 0 and 
U^{t + a', B^>) = 0, P 3 .-a.s. by (3.42). Then subtracting (3.41) from ( 3.51| ) gives 


5< Ex{Ll, - - fi){dy) 

J M 

= / ^x^{a'>{T-t--il}{y))+}{L^ai — L\T-t-'tl)(y))+)^^ ~ 


(3.52) 


We observe that for y E suppz/ one has = 0 (recall c±{t) > a±, t E [0,T)) and then 
for any such y and P 2 ,-a.e. oj E {a' > T — t — 'ip{y)} we have = T — t. Using the 


latter inside (3.52) we get 


/ E3;l{(T'>{T-t-p(y))+} L^iT-t-^ij^yp+^b'idy) <t) 

J M 

which leads to a contradiction. 


(3.53) 

□ 
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Figure 2: Boundary of Rost’s barriers for uniform-like distributions. (Left) = 

0.5 1[o^2](3j)) h = 5 X 10“^ (dashed line) and F'^{x) = 0.5 l[o^o.4](3^) + 0-5 l[o.6,2.2](3;)) h = 10“'^ 
(solid line). (Right) F'^{x) = 0.5 l[o,o.4](3^) + 0-5 1 [ 0 . 6 , 2 . 2 ] ( 2 ^)) h = 4 x 10“^ (solid line), 
h = 2 X 10“^ (dashed line), h = 2 x 10“^ (fine dashed line); as expected the sharpness of 
the jump increases as h fO. 




Figure 3: Boundary of Rost’s barriers: (Left) g = with mean m = 1 and variance 

cj^ = 1 (h = 10“^). (Right) jj, = Exp{X), h = 2 x 10“^, A = 1.5 (solid line), A = 2.5 (fine 
dashed line), A = 3.5 (dashed line). 

4 Numerical solution to the integral equations 

We proved above that s+ and s_ are the unique solutions to the equations ( |2.8[ ). These 
equations cannot be solved analytically but in this section we describe a relatively simple 
algorithm for their numerical solution. 

The convergence of this algorithm has been studied in the literature in several specihc 
examples. For instance in 1201 the authors deal with optimal exercise boundaries for 
American Put and Call options on dividend-paying assets, whereas in Section 5 of [36] 
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the algorithm is applied to inverse hrst passage time problems. A broader overview of 
numerical methods for Volterra equations may be found in [22] as well as in the list 
references of |20] and [36] . 


Now we illustrate the method. As in the previous section, we focus on time reversed 
boundaries b±{t) = s±{T — t) for t G [0, T]. For Z = 0,1,..., N we set ti = Ih with h = T/N. 
Then for fc = 0,1,..., A^ — 1 we consider the following approximations of the left-hand side 
of the integral equations (2.8): 



rb+iu) 

I p{tk,±b±{t),u,y){u - 

-b- (u) 


rb+bi) 

ia){dy)du^h / p{tk,b±{tk),ti,y){iy- 

i=k+r 


p){dy). 


Setting the above equal to zero, with a slight abuse of notation we obtain the approxi¬ 
mating integral equations 


fb+di) 

^ zZ / p{tkMitk)Rhy){T^ - p){dy) =b, (4.1) 

l=k+\ -^-b-dl) 
fb+di) 

^ 2-, / p{tk,-b-{tk),ti,y){iy - p){dy) = 0, (4.2) 

l=k+l d-b-pi) 


which we are now going to solve recursively starting from t]sf = T. 

For k = N—1 and 6±(fv) = b± equations (4.1) and (4.2) become algebraic equations for 
6+(fAr_i) and b-{tN-i), respectively. These can be solved numerically using, e.g., Newton- 
Raphson or bisection method. Next we use the known values of 6+(fAr), 6+(fAr_i), b_{tN) 
and b-{tisf-i) to obtain algebraic equations from (4.1) and (4.2) when k = N — 2. By 
solving the latter numerically we get b^{tj^_ 2 ) and b_{tj^_ 2 ). The procedure is then iter¬ 
ated to compute 6+(tAr), 6+(tAr_i),..., 6+(ti), 6+(to) and b-(tN),b-(tN-i), ■■■,b-(ti),b_(to) 
as approximations of the optimal boundaries 6+ and 6_, respectively, at the points 
T,T — h ,..., h, 0. 


Letting h —)■ 0, arguments similar to those developed in [201132] guarantee convergence 
of the approximations to the actual boundaries. A complete proof of the convergence is 
outside the scopes of this work, however we observe it very clearly in our numerical 
experiments illustrated in Figures 2 and 3. 
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